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MODEL  FOR  NEARFIELD  ELECTROMAGNETIC  SHIELDING 
BY  CYLINDRICAL  SHELLS  OF  COMPOSITE  MATERIALS 

1.  INTRODUCTION 


This  investigation  focuses  on  the  shielding  of  sources  at  low  frequencies  by  cylindrical 
shells  of  composite  materials  characterized  by  complex  permeability  and  permittivity.  In  brief,  the 
magnetostatic  problem  is  solved  by  first  formulating  the  nearfield  source  in  binomial  expansions 
and  the  fields  around  the  cylindrical  shell  in  powers  of  Legendre  polynomials.  1-5  Next,  the 
boundary  conditions  of  the  fields  are  applied  to  the  cylindrical  surface.  The  magnetostatic  solution 
is  then  used  as  input  to  alternate  iterations  of  Faraday's  law  and  Ampere’s  law  to  generate  a 
perturbation  expression  of  the  fields.  The  primary  objective  of  this  process  is  to  develop  a 
mathematical  model  of  the  H,  B,  E,  and  D  fields*  as  they  exist  in  the  presence  of  the  cyhndrical 
shell. 


More  specifically,  the  time-averaged  Poynting  vector  must  be  obtained.  This  requires 
solution  of  Maxwell's  equations  by  iteration,  and  begins  with  calculations  of  the  B  and  H  fields  in 
the  magnetostatic  limit  for  source  models  representing  transformers  and  power  supplies  near  the 
cylindrical  shell.  The  boundary  conditions  for  the  continuity  of  the  tangential  component  of  the  H 
field  and  the  normal  component  of  the  B  field  are  applied  at  the  two  interfaces  of  the  cylindrical 
shell,  as  shown  in  figure  1. 

The  resulting  magnetostatic  fields  are  then  substituted  into  Ampere's  circuital  law  by 
equating  the  curl  of  the  H  field  to  the  current  density  and  the  time  derivative  of  the  displacement 
vector.  The  E  field  can  be  extracted  from  the  current  density  by  using  Ohm's  law.  The  solution  for 
the  E  field  is  generated  from  and,  by  iteration,  is  substituted  into  Faraday’s  law  to  yield  improved 
H  and  B  fields.  The  iteration  proceeds  back  and  forth  between  Ampere's  law  and  Faraday's  law, 
generating  an  infinite  series  of  H  and  B  correction  terms,  as  described  by  Feynman.^  This  series  of 
E  and  H  fields  (truncated  for  the  purpose  of  this  investigation)  will  yield  the  E  and  H  fields  used  in 
the  calculation  of  the  Poynting  vectors.  Because  this  is  a  nearfield  model,  retarded  time  effects  can 
be  ignored. 


*  According  to  Feynman's  definition,  B  represents  the  magnetic  field  vector  and  H  represents  the 
magnetizing  field.  The  term  E  represents  the  electric  field,  and  D  is  the  displacement  field. 
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MEASURING  DEVICE  AT 
POINT  OF  OBSERVATION 


Note:  The  power  supply  field  varies  as  \l  \  f  —  f'  \  and  the  transformer  field  varies  as  1/ 1  r  t  \  ,  where  the  term  r 
is  the  distance  from  the  source  to  the  center  of  the  cylindrical  shell  and  r  is  the  distance  from  the  cylindrical  axis  to 
the  measuring  device  at  the  point  of  observation.  The  coordinates  at  the  point  of  observation  are  always  f  and  ft 


Figure  1.  Representation  of  the  Cylindrical  Shell  (or  Blockage) 
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2.  THE  MATHEMATICAL  MODEL 


As  described  in  section  1,  the  model  is  based  on  an  iteration  process  between  Ampere's  law 


VxH  =  J  +  ^  , 
dt 


where  J  is  the  current  density  and  D  is  the  displacement  field  vector,  and  Faraday's  law 


Vx£=-— , 
dt 


where  V  x  is  the  curl  of  the  E  field  and  t  is  time.  The  H  field  is  related  to  Bhy  B  =  idH,  where 
the  permeability  fx  is  complex  for  a  composite  material.  The  displacement  field  D  is  related  to  E 
by  D  =  £E,  where  E  is  complex  for  a  composite  material. 

Ohm's  law  is  used  to  express  the  current  density  in  terms  of  the  E  field  as 


J  =oE , 


where  cr  is  the  conductivity. 

We  use  the  complex  amplitudes  of  the  fields, 

E  =  ReE^  , 

where  Re  represents  the  real  part  and  E  is  the  complex  amplimde  of  the  electric  field,  to  generate 
the  Poynting  vector 


5=-Re(£xH*)  . 
2 


The  Poynting  vector  is  then  used  to  derive  field  intensities  at  the  point  of  observation: 
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The  iteration  process  can  be  truncated  after  only  a  few  iterations  owing  to  the  presence  of 
low  frequencies.  In  fact,  for  very  low  frequencies,  we  may  only  need  Ampere's  law  to  extract  the 
E  field  from  the  H  field  (which  has  been  derived  from  the  magnetostatic  model). 
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3.  THE  POWER  SUPPLY  SOURCE  IN  THE  PRESENCE  OF  A 
CYLINDRICAL  SHELL  OF  COMPOSITE  MATERIAL 
IN  THE  MAGNETOSTATIC  LIMIT 


THE  MATHEMATICAL  FORM  FOR  THE  l/|f-r'|  SOURCE 

For  the  power  supply  source,  the H fields  are  given  in  three  regions:  r>b,  a<r<b,  and 
r<a,  where  a  and  b  are  the  inner  and  outer  radii  of  the  cylindrical  shell,  respectively. 


For  r  >b,  the  H  field  is  given  by 


H  = 


r-r\ 


\r-r 


£=0 


(1) 


where  Pg{cos6)  are  the  Legendre  polynomials. 


For  a<r<b. 


(2) 


and  for  r<a. 


2 

// =  -V^5^/p^[cos(0)]. 

^=0 


(3) 


THE  l/|r-ri  SOURCE  IN  COMPONENT  FORM  AND 
ITS  BINOMIAL  EXPANSION 


The  B  field  for  the  l/|f  -  f'\  can  be  expressed  as 


B,  (r-n 

\f -f'\  If-f'l  ’ 
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where  the  second  factor  is  the  unit  vector  extending  from  the  source  to  the  observation  point. 
Taking  the  dot  product  of  B  with  the  unit  vectors  in  the  radial  and  angular  directions  of  the  physical 
blockage  cylindrical  coordinate  system  yields 


B  = 


Bo  ir-r) 


{r-r'cos6) 


(4) 


and 


R  -  B,  f  r'\*n-n 


(5) 


where  B^  is  the  source  strength. 


Expanding  1  / 1  f  -  f'  1^  in  the  binomial  expansion  yields 


r’ 


r'^  2r'  ^ 

1 - r-H - COS0 

V  r"  r  y 


for  f>f' 


and 


r  J2- 


\r-ff  r'^ 


K  2r 

1 - ^+— COS0 

r  ^  r 


for  f  <f' . 


(6) 


(7) 


THE  MATHEMATICAL  SOLUTION  OF  THE  B  AND  H  FIELDS 
FOR  THE  l/|r-r'|  SOURCE 


We  substitute  equations  (6)  and  (7)  into  equations  (4)  and  (5)  and  the  result  into  the  source 
terms  of  equations  (1),  (2),  and  (3).  By  applying  the  boundary  conditions  for  the  B  and  H  fields 
to  the  resulting  equations  and  setting  the  coefficients  of  the  powers  of  the  sines  and  cosines  equal 
to  zero,  we  generate  10  equations  and  10  unknowns,  which  are  solved  to  obtain  the  coefficients  of 
expansion  with  the  MathematiccP  software  (see  appendix). 

The  boundary  atr  =  b  has  two  cases:  f  <f'  and  f  >  r'.  For  the  continuity  of  Br ,  the 
coefficients  of  (cos  0)®  are 


A. 


^-72 

\  ^  J 


(8) 


6 


and 


bJ. 

b - 

V  b  j 


13  1  1  3  y 


the  coefficients  of  (cos  6^  are 


B, 


^3h^ 


-r 


=  -2ai-^-2/xyi-^  for  f<r' 
b  b 


(10) 


and 


b‘ 


/„/3  ^ 

-  +  r' 
\o  j 


b  b 


(11) 


and  the  coefficients  of  (cos0)^  are 


-IbB^  9  1  ooL^y,- 

-75^  =  :7«2  — +  for  r<r 


2 


2"  h" 


(12) 


and 


^^^3^  =  f«2-T  +  3AtB2^-|/x^  for  f>f' 
b  Z  r  Z  o 


(13) 


For  the  continuity  of  Hq  at  r  =  b,  there  are  no  terms  in  the  coefficient  for  (cos  0)®.  For  the 
coefficient  of  (sin)i, 


B. 


(14) 


and 


r  r'^\ 

r'-  — 

K  J 


-^--R  -II 
b^  b^ 


for  f>f' 


(15) 


For  the  coefficient  of  (sin  Bcos  0)^, 


IBob  3a,  „o.2  3y,  . 

-^  +  -^  =  3l5jy^  +  -^  for  r<r' 


(16) 


7 


and 


3a,  oo  1.2  372  r 
-^+-p^  =  3/?2^  +-^  for  r, 


(17) 


For  the  continuity  of  at  r  =  a,  the  coefficient  of  (cos0)®  is 


a 


In' 


(18) 


the  coefficient  of  (cos0)^  is 


m(^A- 


^1  =  5- 


(19) 


2 

and  the  coefficient  of  (cos0)  is 


(20) 


For  the  continuity  of  H  at  r  =  a,  the  coefficient  of  (sin)i  is 


A^+ri^=A«  ’ 


(21) 


and  the  coefficient  of  (sin  6  cos  0  )i  is 


A«^  +  r2^  =  • 


(22) 


The  radii  a  and  b  are  fixed  known  numbers,  and  the  unknowns  are 


Uq,  Up  <X2,Pi,p2’  7o>  7i.  72'  <5;i,and  ^2  . 


with  the  following  solutions: 


8 


+  +  18/i<3V  +  +  6juV) 


6X2X4/ 


-6b^X2X2  - S^iXzXsr' -  6ii^X2Xs/  - 6bX2X2/^  +  6b'^X4X^/^] 


6X2X4/ 


B,{-3b^  +  /^) 


2/^ 


*  2X2/^  ’ 


2ByX, 

3X,/^ 


10  B^lLb^X, 

^  3X2  r'^  ’ 


which  were  solved  numerically  using  the  Mathematical  software  (see  appendix). 

Substituting  the  above  Legendre  polynomial  coefficients  back  into  equations  (1),  (2),  and 
(3)  yields  foTr>b 


H  =  r  -2 - ^ - (r-r'cos0)--^cos0+-^  +  -^^(3cos^  0-l] 

r+/^-2r/  r  r  2r  ^  \ 


Bnr'sinS  a. 


- ^sm0  +  — x(6cos0sin0) 


9  9 - A 

r  +/  -2r/  r  2r 


for  b>r>a 


^  =  n  n 7  " ^1 ^ 7^  +  [/i 7  +  ^P^rpcos^  d-l) 


67,  . 

_ ^  01  r 


+  0l)3j+-f  sin0-)52r  +  — 52-sin0cos0  ; 


and  for  r  < 


H  =  r  ^^cos0-'^%(3cos^0-l]  +  ^  +0f'^sin0  +  ^^sin0cos0 
[^3  4^4  \  I  ^2J  / 


These  equations  have  complex  permeability. 
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4.  THE  TRANSFORMER  SOURCE  IN  THE  PRESENCE  OF  A 
CYLINDRICAL  SHELL  OF  COMPOSITE  MATERIAL 
IN  THE  MAGNETOSTATIC  LIMIT 


THE  MATHEMATICAL  FORM  FOR  THE 
SOURCE 

For  the  transformer  source,  the  H  fields  are  given  in  three  regions:  r>b,  a  <r<b,  and 
r<a,  where  a  and  b  are  the  inner  and  outer  radii  of  the  cylindrical  shell,  respectively.  For  r>b, 
the  H  field  is  given  by 


r-r  I  |r-r  | 


for  a<r<b. 


and  for  r<a. 


2 

H  =  [cos(e)] . 

.^=0 

Expanding  equations  (23),  (24),  and  (25)  yields  for  r  >  h 


(23) 


(24) 


(25) 


H  =  r 


Bo 


+  r'^  -  2rr'J 


-(r  -  r'cos0)  -  ^cos0  + 


CXq 


3a2 

2r^ 


(3cos^0-l) 


+  0 


By^'sin0 


+  r'^  -  2rr'j 


— — sin  0  + (6  cos  0  sin  0) 


2r 


ioT  b>r> a 


11 


H  =  r 


7i  ^  -  a]cos0  +  ^  +  (^72  ^  +  2/?2r)(3cos2  0-l) 


+  0 


Bi  +  ^\in0  -  J32r +^^sin0cos9 
r  J  r 


and  for  r  <  ^? 


H  =  f\ 


2Si  33^  2 

—!^cos6 - 2-(3cos^  0-1)+ 

r  4r 


A 

+  6 


’^1  .  . 
^sine  + 

r 


sin  0  cos  6 


THE  l/|r-r'f  SOURCE  IN  COMPONENT  FORM  AND 
ITS  BINOMIAL  EXPANSION 

The  B  field  for  the  l/|f  -  f'\  source  can  be  expressed  as 


where  the  second  factor  is  the  unit  vector  extending  from  the  source  to  the  observation  point. 
Taking  the  dot  product  of  B  with  the  unit  vectors  in  the  radial  and  angular  directions  of  the  physical 
blockage  cylindrical  coordinate  system  yields 


B  = — - ^(r-r'cos0) 


(26) 


and 


(27) 


\r-r' 


r  —  r 


Expanding  l/|f  -  f'f  in  the  binomial  expansion  yields 


/  «-l- 

1—  — /|4  4 

1-2 

^-^cos^l 

|r  -  r  1  r 

for  r>f' 


(28) 
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and 


THE  MATHEMATICAL  SOLUTION  OF  THE  B  AND  H  FIELDS 
FOR  THE  l/|r-r'p  SOURCE 

We  substitute  equations  (28)  and  (29)  into  equations  (26)  and  (27)  and  the  result  into  the 
source  terms  of  equations  (23),  (24),  and  (25).  By  applying  the  boundary  conditions  for  the  B  and 
H  fields  to  the  resulting  equations  and  setting  the  coefficients  of  the  powers  of  the  sines  and 
cosines  equal  to  zero,  we  generate  the  10  equations  and  10  unknowns,  which  are  solved  to  obtain 
the  coefficients  of  expansion.  These  coefficients  are  then  substituted  back  into  equations  (1),  (2), 
and  (3)  to  yield  a  numerical  solution. 

The  boundary  at  r  =  b  has  two  cases:  f  <  T  and  r  >  f'.  For  the  continuity  of  4 ,  the 
coefficients  of(cos0)°  are 


B. 


1- 


2Kb 


2.  - 


for  r<r' 


(30) 


and 


b- 


^  1  2 
V  ^  y 


cCq  3  f  0^2^  _ 


— r+- 
b^  2 


^4  i  =  -/^ro 


(-L] 

W) 


for  f>r\  (31) 


the  coefficients  of  (cos  0)^  are 


B, 


Ub^  -  2b 
— — b+  — 


2V 


I  for  f<r 


(32) 


and 


B, 


f  r'^\  2a  (' 

3A-'  +  Tr|--^  =  M-^7if^l  for  r>r'; 


J 


(33) 


and  the  coefficients  of  (cos  d'f  are 


13 


for  f<f' 


(34) 


and 


9«2  ,3/^ 
2 


for  r>r'  . 


For  the  continuity  of  at  r  =  b,  there  are  no  terms  in  the  coefficient  for  (cos  0)°. 
coefficients  of  (sin)^ 


A 

r' 


r'- 


CL, 


1  _ 


Pib  +  Yi^l  for  r<r' 


and 


B, 


r  — ■ 


2r 


'3\ 


a, 


1  _ 


A&+ri4^ 

b  J 


for  f>f' 


For  the  coefficients  of  (sin  6  cos  6)^, 


For  the  continuity  of  B^dX  r  —  a,  the  coefficient  of  (cos  6)^  is 


the  coefficient  of  (cos  6^  is 

A-ri^=A; 

and  the  coefficient  of  (cos  6)^  is 


(35) 

For  the 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 
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(42) 


f(2Aa-n^) 


=  35^a  . 


For  the  continuity  of  at  r  =  a,  the  coefficient  of  (sin)^  is 

V  a  J 

and  the  coefficient  of  (sin  6  cos  Oy  is 


(43) 


(44) 


Solving  the  10  equations  and  10  unknowns  with  p.  complex,  which  represents  a  composite 
material,  yields  the  magnetostatic  solution  (calculated  with  Mathematical  as  shown  in  the 
appendix). 
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5.  CONCLUSIONS 


This  research  develops  a  new  technique  to  model  fields  from  transformers  and  power 
supplies  at  low  frequencies  in  the  nearfield  region  shielded  by  cyhndrical  shells  of  composite 
materials.  The  first  step  in  the  process  formulates  the  magnetostatic  problem  by  solving  Laplace's 
equation  in  terms  of  binomial  and  Legendre  polynomial  expansions.  With  the  solution  to  the 
magnetostatic  problem  used  as  a  basis.  Ampere's  law  generates  an  E  field  with  complex 
permittivity.  This  E  field  becomes  the  basis  for  a  new  H  field  in  Faraday's  law.  A  perturbation 
expansion  is  then  generated  for  the  B,  H,  E,  and  D  fields.  The  boundary  conditions  are  applied  to 
these  new  fields,  and  the  coefficients  of  expansion  are  calculated.  The  intensities  and  Poynting 
vectors  are  also  obtained  in  the  various  regions.  Plots  of  these  intensities  versus  the  real  and 
imaginary  parts  of  the  permeability  and  the  real  and  imaginary  parts  of  the  permittivity  will  be 
provided  in  a  future  report. 


17/18 
Reverse  Blank 


6.  REFERENCES 


1 .  G.  Arfken,  Mathematical  Methods  for  Physicist,  Third  Edition,  Academic  Press,  New  York, 
1985,  pp.  637-707. 

2.  L.  Bailin  and  S.  Silver,  "Exterior  Electromagnetic  Boundary  Value  Problems  for  Spheres  and 
Cones,"  IEEE  Transactions  on  Antennas  and  Propagation,  vol.  AP-4,  no.  1, 1956,  p.  15. 

3.  D.  Jackson,  Classical  Electrodynamics,  First  Edition,  John  Wiley  &  Sons,  Inc.,  New  York, 
1962. 

4.  T.  Anderson,  "Magnetic  Field  Intensity  Due  to  a  Line  of  Current  of  Arbitrary  Orientation  in 
Cartesian  Coordinates,"  NUWC-NPT  Technical  Memorandum  941145,  Naval  Undersea 
Warfare  Center  Detachment,  New  London,  CT,  23  November  1994. 

5 .  J.  Reitz  and  F.  Milford,  Foundations  of  Electromagnetic  Theory,  Addison-Wesley  Publishing 
Company,  Inc.,  Reading,  MA,  1967. 

6.  R.  P.  Feynman,  Lectures  in  Physics:  Volume  II,  Addison-Wesley  Publishing  Company,  Inc., 
Redwood  City,  CA,  1962. 

7 .  S.  Wolfram,  Mathematica:  A  System  for  Doing  Mathematics  by  Computer,  Second  Edition, 
Addison-Wesley  Publishing  Company,  Inc.,  Redwood  City,  CA,  1991. 


19/20 
Reverse  Blank 


APPENDIX 


THE  USE  OF  MATHEMATICA  TO  OBTAIN  A  SPECIFIC  SOLUTION 
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Reverse  Blank 


■  Description 

This  Mathematica  program  solves  the  magnetostatic  problem  of  a  power  supply  interacting  in  the  nearfield  with 
cylindrical  shell  of  composite  material  that  has  complex  permeability. 
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Set  all  the  coefficients  equal  to  0 
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Simplify  the  final  equations  for  the  coefficients  of  the  Legendre  Polynomials 

Simplify [Coeffs] 
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Substitute  the  Coefficients 
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